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ABSTRACT 

The 229 GHz (A1.3mm) radio emission from Orion-KL was mapped with up to 0.14" angular 
resolution with CARMA, allowing measurements of the flux densities of Source I ('SrcF) and the 
Becklin-Neugebauer Object (BN), the 2 most massive stars in this region. We find integrated flux 
densities of 310 ± 45 mJy for Srcl and 240 ± 35 mJy for BN. Srcl is optically thick even at 229 GHz. 
No trace of the H30a recombination line is seen in its spectrum, although the i>2=l, 5(5,0)-6(4,3) 
transition of H2O, 3450 K above the ground state, is prominent. Srcl is elongated at position angle 
140°, as in 43 GHz images. These results are most easily reconciled with models in whic h the radio 
emissi on from Srcl arises via the H free-free opacity in a T < 4500 K disk, as considered by lReid et al.1 
(|2007| ). By contrast, the radio spectrum of BN is consistent with p + /c~ free-free emission from a dense 
(n e ~ 5 x 10 7 cm -3 ), but otherwise conventional, hypercompact HII region. The source is becoming 
optically thin at 229 GHz, and the H30a recombination line, at Vlsr=23.2±0.5 km s _1 , is prominent 
in its spectrum. A Lyman continuum flux of 5 x 10 45 photons s _1 , consistent with that expected from 
a B star, is required to maintain the ionization. Supplementary 90 GHz observations were made to 
measure the H41a and H42a recombination lines toward BN. Published 43 and 86 GHz data suggest 
that Srcl brightened with respect to BN over the 15 year period from 1994 to 2009. 
Subject headings: ISM: individual(Orion-KL) — radio continuum: stars — radio lines: stars — stars: 
formation — stars: individual (Bccklin-Neugebauer-Objcct) 



1. INTRODUCTION 

The Kleinmann-Low Nebula in Orion is well-known as 
the n earest region of high mass star form ation, 415 pc 
away (jMenten et all 120071 : iKim et all 120081 ). It contains 
at least 2 massive young stars. One of these, the Becklin- 
Neugebauer Object (B N) , has been studied extensively at 
infrared wavelengths dScoville et al.|[l983l ); it is thought 
to be a B star. The other object, Source I (hereafter, 
'SrcF), is so heavily obscured by foreground dust that it 
is not directly visible in the infrared, although light re- 
flected by the sur rounding nebulosity provides a glimps e 
of its spectrum ([Morino et alj|1998t iTesti et al.M201Ql) . 
Srcl is noteworthy because it is surrounded by a cluster 
of SiO masers, one of the few cases in which SiO masers 
are associated with a young star. It also is known to 
drive a bipolar outflow into the surrounding molecular 
cloud (jPlambeck et alJl2009f ). 

Remarkably, proper motion measurements show that 
Srcl and B N are recoiling from one another at 35- 
40 km s" 1 (iRodriguez et al.l [20051 : iGomez et all [200l 
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Goddi et al. 2011). Tracing the motions backward, 
Goddi et al.l (| 2 1 lh find that 560 years ago the projected 



separation of these two stars in the plane of the sky 
was just 50 ± 100 AU. An extensive system of 'bullets,' 
bow sho cks, and 'fingers ' vis i ble in lines of H2, Fell, 
and CO (jAllen fc Burtonlfl99lHZapata et aJJ[201lt >. also 
is centered on Orion-KL. Proper motion measurements 
suggest that the fingers w ere created by a n explosive 
event 500-1000 years ago (jDoi et al.l 120021 IBallv et al.l 
1201 ll) . Thus, the curren t ly favored paradigm for Orion- 
KL jGomez eTaTl [20081 IBallv et al.l \20TVt IGoddi etall 
120111) postulates that Srcl and BN were ejected from a 
multiple system approximately 500 years ago, and that 
the ejection of the stars unbound the surrounding gas 
and circumstcllar disks, creating the fin ger system. 

BN 's luminosity of (0.8-2.1) x 10 4 L Q (|De Buizer et al.l 
I2012D suggests tha t it is a 10 to 15 M Q star 
( Schaller et al.lll992[ ). Estimates of the mass of Srcl are 
conflicting. Srcl is recoiling at about half the speed 
of BN in the rest fr ame of the Orion Nebula Clus- 
ter (jGoddi et al.ll201lT ) ; conservation of momentum then 
suggests that its mass is roughly 20 M Q . On the other 
hand, if the SiO masers near Srcl are in Keplerian ro- 
tation about the star, the in ferred central mass is only 
7 M Q ([Matthews et al.lfeOlOl ): this can be interpreted as 
a lower limit if the SiO-cmitting gas is supported in part 
by radiation or magnetic pressure. Srcl could well be 
a compact binary with a scmimajor axis of < 10 AU 
(|G6mez et al.|[2008[ ). 

At cm and mm wavelengths the continuum spectra of 
BN and Srcl are believed to be dominated by free-free 
emission from ionized circumstcllar gas that is at least 
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parti ally optically thick up to 100 GHz (Plambcc k et al.l 
1 1995D . Observations at higher frequencies, where this 
emission should become optically thin, can be used to es- 
timate the Lyman continuum fluxes from the stars, bet- 
ter constraining their masses. Such observations must 
be made with subarcsecond angular resolution in order 
to distinguish the compact circumstellar emission from 
bright but extended dust and molecular line emission 
from the surrounding molecular cloud. 

Here we describe high angular resolution 229 GHz ob- 
servations of Orion-KL with the Combined Array for Re- 
search in Millimeter- Wave Astronomy (CARMA) that 
cleanly resolve the emission from both stars. Surpris- 
ingly, we find that the emission from Srcl is optically 
thick even at 229 GHz, and we fail to find the H30a re- 
combination line in its spectrum. We argue that these 
results favor models in which the radio emission from 
Srcl originates from the H~ opacity in a T < 4500 K 
disk. 

2. OBSERVATIONS 

The spectra of Srcl and BN presented here were ob- 
tained from a single night's observation of Orion made 
with the CARMA B-array on 2011 January 06. Four- 
teen antennas were used; projected antenna separations 
ranged from 50 to 700 kA. The phase center was 
05 h 35 m 14?505, -05°22'30'/45 (J2000). The weather was 
excellent - the atmospheric opacity at 225 GHz was 0.10, 
and the rms atmospheric phase fluctuations were 65 fim 
on a 100-m baseline. Double sideband system tempera- 
tures for the dual-polarization receivers were 100-200 K, 
scaled to outside the atmosphere. For each polarization 
the correlator was configured to observe three overlap- 
ping 500 MHz windows (covering the frequency range 
225.8-227.2 GHz in the receivers' lower sideband and 
231.4-232.8 GHz in the upper sideband, with 12.5 MHz 
wide channels), plus a single 250 MHz window (covering 
226.65-226.90 GHz in the lower sideband and 231.75- 
232.00 GHz in the upper sideband, with 3.125 MHz res- 
olution) . Upper and lower sideband signals are separated 
in the correlator by phase-switching. 

The antenna gains were derived from observations of 
the calibrators 0423-013 (2.4 Jy) and 0607-085 (1.05 Jy) 
that were interleaved with the Orion scans every 10 min- 
utes. To minimize decorrelation from atmospheric phase 
fluctuations, a self-calibration interval of 30 seconds was 
used. The calibrator flux densities were established from 
observations of the primary flux standard, Uranus. Be- 
cause Uranus is heavily resolved on these baselines, only 
antenna pairs with separations < 200 kA were used for 
flux calibration. Antenna gain solutions were applied 
to the Orion data after smoothing to a 10 min interval 
(scalar averaging the amplitude corrections). Only data 
with antenna separations greater than 250 kA were used 
for the Orion maps, yielding a 0.28" x 0.21" synthesized 
beam at PA -27°. 

With 0.25" angular resolution, almost all spectral line 
emission from the Orion Hot Core and surrounding 
molecular cloud is resolved out. Nevertheless, to ensure 
that continuum flux densities were not contaminated by 
residual line emission or absorption, we made maps of 
all the spectral channels, and flagged those (28 out of 
234) with anomalously high (> 2x average) rms noise; 
the noisy channels coincided with strong spectral lines. 



We generated a primary-beam corrected continuum map 
from the remaining channels, using multifrcquency syn- 
thesis to avoid radial smearing. 

To help correct the Orion data for atmospheric decorre- 
lation, a single phase-only self-calibration was performed 
on the source itself. This means that after cleaning the 
continuum map in the normal way, the brightest pixels 
in the clean component list were used as a source model 
to refine the antenna-based phases. A time interval of 2 
minutes was used, short enough to track the most signif- 
icant atmospheric phase variations, but long enough to 
obtain a signal to noise ratio of >1 on each of the 91 avail- 
able baselines, allowing a robust least squares fit for the 
13 antenna-based phases. Applying the self-calibrated 
phase corrections increased the source flux densities by 
about 50%. The measured rms noise in the final 0.25" 
resolution map is 4 mJy beam -1 . 

Even higher resolution images of Orion were ob- 
tained in 2009 February using the CARMA A-array, 
with antenna separations of up to 1.7 km (1350 kA). 
These observations utilized older, single-polarization re- 
ceivers and a 1.5 GHz bandwidth correlator. For these 
data the CA RMA Paired Ant enna Calibration System 
("C-PACS," iPerez et al.l [20l0| ) was used as an adap- 
tive optics scheme to correct for blurring by atmospheric 
phase fluctuations. C-PACS placed the 8 CARMA 3.5- 
mctcr telescopes next to a subset of the 10-m and 6- 
m telescopes, including those at the extremes of the ar- 
ray. As the 6-m and 10-m telescopes observed Orion at 
229 GHz, the 3.5-m antennas observed the nearby (1.6° 
away) quasar 0541-056 at 31 GHz in order to monitor the 
atmospheric phase variations above each station. The 
atmospheric phase delay to each 3.5-m telescope was de- 
rived by self-calibration on 0541-056 every 12 sec. These 
antenna-based phases were scaled up by the ratio of the 
observing frequencies (229/31 = 7.4) and applied to the 
mm data. The C-PACS calibration accuracy is limited 
by the angular separation of the cm calibrator and the 
mm source, which causes the cm and mm beams to probe 
different paths through the atmosphere; by the signal to 
noise ratio on the 31 GHz calibrator; and by the airmass. 

A continuum map was generated from the A-array data 
using only baselines longer than 250 kA. The synthesized 
beam was 0.15" x 0.13" at PA 14°. Applying the C-PACS 
corrections more than doubled the signal to noise ratio of 
this map. As with the B-array map, we then performed a 
single phase-only self-calibration, with a 2 minute inter- 
val, to help correct for residual atmospheric phase fluc- 
tuations; this doubled the flux densities of the compact 
sources. The measured rms noise in the final 0.14" reso- 
lution map is 2.6 mJy beam -1 . 

The integrated flux densities measured from the 0.14" 
and 0.25" resolution maps were 240 and 255 mJy for BN, 
310 and 370 mJy for Srcl. Although the discrepancy in 
the two Srcl flux densities could conceivably be caused 
by time variability of the source (cf. section I4.2.3P , it 
is more likely attributable to differences in the sampling 
of visibilities by the A and B arrays, which lead to dif- 
ferent sidelobe structure from the adjacent Orion Hot 
Core. We will use the A-array flux densities for the re- 
mainder of this paper because the higher resolution data 
should more effectively filter out emission from the Hot 
Core. We estimate that the absolute flux scale is accu- 
rate within ±15%. 
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Figure 1. (left) 229 GHz continuum maps of Orion-KL from CARMA C- and B-array data. Black contours show the 0.83" resolution 
C-array map; the contour levels are —100,-50, 50, 100, 200, 300, and 400 mjy/beam. Negative contours are attributable to poorly 
sampled extended emission. Red contours show the 0.25" resolution B-array image; contour levels are 20, 40, 60, 80, 100, 150, 200, and 
250 mjy/beam. Extended emission is almost entirely filtered out in this map; the most negative value is —30 mjy/beam. (right) Spectra 
of BN and Srcl near 232 GHz obtained from the 0.25" resolution image. Dashed lines indicate the continuum levels. The spectra have 
been Hanning smoothed to 30 km s — 1 velocity resolution. The H30a line is detected toward BN, the H2O r>2=l, 5(5,0)-6(4,3) line toward 
Srcl. Arrows in the Srcl spectrum indicate weak absorption features from transitions of ethyl cyanide, dimethyl ether, and methanol. 



Finally, supplementary 3mm observations were made 
with the A-array in 2009 February to obtain spectra 
of the RAla (92.034 GHz) and H42a (85.688 GHz) re- 
combination lines toward BN. Both of the recombination 
lines and the strong SiO maser at 86.243 GHz were mea- 
sured simultaneously. The data were self-calibrated on 
the maser using a 10 sec interval. The spectra of the two 
recombination lines were averaged together to improve 
the signal to noise ratio. 

3. RESULTS 

The left hand panel in Figure Q] shows the 0.25" res- 
olution continuum map from the B-array observations, 
overlaid on a 0.83 " resolution image obtained with the 
CARMA C-array (jEisner et all 12008]). Srcl and BN are 
clearly distinguishable only in the higher resolution map; 
extended emission from dust and molecular lines domi- 
nate the lower resolution image. 

The upper sideband spectra of BN and Srcl are shown 




Figure 2. Images of BN and Srcl at 229 GHz from the CARMA 
A-array data. Each box is 0.6" on a side. The lowest contour and 
the contour interval is 15 mjy beam -1 . The rms noise level is 2.6 
mjy beam -1 . The 0.15" X 0.13" synthesized beam is shown by 
a hatched ellipse in the upper left corner of eac h image. Srcl is 
extended at PA 14 0° , as in 43 GHz VLA images <Reid et al.|[2007l ; 
IGoddi et afll20TlD . 



in the right hand panels of Figure [TJ The H30a recombi- 
nation line at 231.901 GHz is prominent in the spectrum 
of BN, but is not detected toward Srcl. The bright line at 
232.687 GHz in Srcl's spectrum is the u 2 =l, 5(5,0)-6(4,3) 
transition of H2O. Weak absorption features toward Srcl 
are due to transitions of ethyl cyanide, dimethyl ether, 
and methanol. No emission lines were detected toward 
cither source in the lower sideband. 

Figure [5] displays the 0.14" resolution images of BN 
and Srcl from the A-array. BN is nearly circularly sym- 
metric, while Srcl is elongated. Table Q] summarizes the 
results of Gaussian fits to the positions, flux densities, 
and deconvolved sizes. The 229 GHz source positions 
match the 43 GHz p ositions measured in 2009 January by 
IGoddi et all (|2011l) within 0.015". The deconvolved size 
of Srcl at 229 GHz, 0.20" x 0.03" at PA 140°, is similar 
to th e sizes measured at 8.4 GHz (0.19"x < 0.15" at PA 
136°: iGomez et alJl200l and at 43 GHz (0.23" x 0.12" 
at PA 142°: IGoddi et al.lfeOllh . although the source ap- 
pears to be thinner along its minor axis at 229 GHz. 

The flux density we measure for BN in the 0.14" reso- 
lution map i s roughly 2x higher than th e value recently 
reported by iGalvan-Madrid et al.l ([2012D from an anal- 
ysis of ALMA Band 6 Science Verification data. These 
early ALMA data had a 1.3" x 0.6" synthesized beam. 
We suspect that the flux density of BN is depressed in the 
ALMA map because of negative sidelobes from poorly 
sampled extended emission. The same problem affects 
the 0.83" resolution CARMA image shown in Figure [TJ 
the peak intensity at the position of BN in this low reso- 
lution image is ~ 125 mjy/beam. The higher resolution 
CARMA maps are free of deep negative sidelobes, allow- 
ing more reliable flux density measurements of compact 
sources. 
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Table 1 

Source parameters for Srcl and BN at 229 GHz 



source 


RA 


DEC 


peak flux density 


integrated flux density 


deconvolved size 




( h m s ) 


(° ' ") 


(mJy) 


(mJy) 




Srcl 


05 35 14.514 


-05 22 30.59 


170 ± 25 


310 ± 45 


0.20" x 0.03" at PA 140° 


BN 


05 35 14.109 


-05 22 22.73 


215 ± 30 


240 ± 35 


0.06" x 0.04" at PA 90° 



Note. — All parameters were measured from the 0.14" resolution map. Positions are for epoch 2009.1. 




freq (GHz) freq (GHz) 

Figure 3. Radio spectra of BN and Srcl. Squares indicate continuum flux densities; circles, flux densities at the peak of recombination 
lines. CARMA results presented in this paper are shown as red in the online edition. Solid and dashed lines show, respectively, the 
continuum and recombination line flux densities predicted by the models described in sections 4.1.1 and 4.2.1. 



Figure [3] shows the radio continuum spectra of BN and 
Srcl from cm to submillimeter wavelengths. Numerical 
values are given in Table [5J For BN we also plot the peak 
intensities of the H30a and H41/42a recombination lines 
measured with CARMA, and the H 53a line measured 
with the VLA ([Rodriguez et al.ll2009D . 

4.1. BN 

The continuum spectrum of BN is typical of free-free 
emission from a hypercompact HII region. The flux den- 
sity scales as v , indicating that optically thick emission 
has a larger angular extent at low er frequencies, either 
becau se the ionized gas is clumpy ([Ignace fe Churchwelll 
|2004[) or bec ause the electron densit y declines smoothly 
with radius (|Wright fc Barlowlll975[ ). The spectrum be- 
gins to flatten at 229 GHz, indicating that the source is 
becoming optically thin at higher frequencies. 

4.1.1. Model spectra 

To model the spectrum, we assume that the emission 
region is spherically symmetric, with a uniform electron 
temperature T e and an electron density that is constant, 
n e (r) = hq, inside core radius r c , then declines as n e (r) — 
no (r/r c ) a at larger radii. For each projected radius x we 
numerically integrate along the line of sight to compute 
the emission measure EM(x) = frv^dz pc cm~ 6 . The 



brightness temperatures at this projected radius in the 
continuum and at the peak of a recombination line are 
then T c {x) = T e [l - exp(-r c )] and T L {x) = T e [l - 
cxp(— tc — tl)], wher e the continuum and li ne opacities 
are approximated as (Rohlfs & Wilson 2000) 

t c (x) = 8.23 x 10- 2 T- 1 - 36 ^ghz EM{x) (1) 
t l (x) = 1.92 x 10 3 T- 2 - 5 EM(x) A^ kI j z . (2) 

We integrate the brightness temperatures over circular 
annuli to obtain flux densities, assuming a distance to 
the source of 415 pc. 

The model BN spectra shown by smooth curves in 
Figure |3] assumed T e = 8000 K, n Q = 5 x 10 7 cm" 3 , 
r c = 7.4 AU, a = —3.5, and recombination line veloc- 
ity widths of 30 km s _1 for all transitions. The model 
does not take into account systematic velocity gradients 
or pressure broadening of the recombination lines. The 
model parameters were determined by trial and error, 
and are not a formal fit to the data. To first order, how- 
ever, the turnover frequency sets uq, the spectral slope 
sets a, and the absolute flux densities set r c . 

4.1.2. Excitation parameter 

The model parameters above yield an excitation pa- 
rameter U = r c (n e Hjj) 1//3 = 5 pc cm~ 2 , which serves as 
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Table 2 

Integrated flux densities for Srcl and BN 



frequency 


S(SrcI) 


S(BN) 


epoch 


References 


(GHz) 


(mjy) 


(mjy) 






4.8 




2.2 ± 0.8 


1990 


1 


8. 1 


1.1 ± 0.2 


4.5 ± 0.7 


1994.3 


2 


8.4 


1.2 ± 0.1 


4.8 ± 0.1 


2006.4 


3 


15 




9.1 ± 0.6 


1981.6 


i 


15 




9.4 ± 1.0 


1983.7 


■1 


15 


1.54 ± 0.18 


6.19 ± 0.19 


1986.3 


5 


15 


1.6 ± 0.4 


6.8 ± 1.5 


1990 


1 


22 




13.6 ± 1.1 


1983.7 


i 


22 


5.7 ±0.9 


16.57 ±0.7 


1991.5 


6 


13 


13 ± 2 


31 ±4.7 


1994.3 


2 


43 


10.8 ±0.6 


28.0 ±0.6 


1994.9 


7 


43 


13 


26.4 ±0.7 


2000.9 


3.8 


43 


14.5 ±0.7 


28.6 ±0.6 


2007.9 


9 


43 


11 ± 2 


23 ±2 


2009.0 


10 


86 


34 ±5 


84 ± 10 


1995.0 


11 


89 


50 ±5 


91 ±9 


2009.0 


12 


229 


310 ± 45 


240 ± 35 


2009.1 


13 


348 


320 ± 48 




2004.1 


11 


690 


6700 ± 3200 




2005.1 


15 



Re ference s. — (11 IFelli et al.l I1993bl ; (21 IMenten fc Reidl 
1995; (31 IGomez et alJ 1 20081 ( 41 IGarav et al.l 1987; 

(5) IFelli et alJ I1993al ; (6) IForbrich et"alT 2008; 

(7) IChandler L Wood! [19971 ; (8) IReid etaLl 2007; 
(9) Rodriguez et al. 2009; (10) Goddi et al. 2011; 
(111 Plambeck et al] [19951 ; (12) IFriedel fc W idicus Weaver 
120111 (assuming 10% absolute c alibr ation accuracy); (13 ) this 
paper; (14l lBeuther et al.ll200l (151 IBeuther et aljflioBI 

a measure of the flux of ionizing photons from the cen- 
tral star. This corresponds to the excitation parameter 
expected for a main sequence B star with Lyman contin- 
uum flux 5 x 10 45 photons s -1 and total luminosity (0.5 — 
1) x 10 4 L Q (IPanagialll97l . The ionizing flux is c ompa- 
rable to the value derived by iScoville et alJ (|1983| ) from 
infrared data, and the luminosity is c omparable with 
the bo lometric luminosity obtained by iDe Buizer et al.1 
(|2012f) from 6-37 /im SOFIA observations. 

4.1.3. Recombination lines 




-100 -50 50 100 



V LSR (km/sec) 

Figure 4. Hydrogen recombination lines observed toward BN. 
H30« and H41/42« (the average of th e two spectra) were ob served 
with CARMA; H53a, with the VLA IjRodrfguez et al.ll2009T ). The 
H53a flux densities are scaled up by 4 for better visibility. Smooth 
curves show Gaussian fits to the spectra. 



Figure 2] displays the recombination line profiles to- 
ward BN. H30a and H41/42a (the average of the H41a 
and H42a lines) were observed wi th CARMA, H53a with 
the VLA (jRodriguez et al.ll2009D . The H30a spectrum 
was measured using the 250 MHz wide correlator win- 
dow with 4 km s _1 channel spacing; it has been Hanning 
smoothed to 8 km s _1 resolution. Line parameters and 
uncertainties derived from Gaussian fits to the spectra 
are summarized in Table [3l The H53a line is at slightly 
lower velocity and is noticeably wider than the higher 
frequency lines. 

For the electron densities n e > 10 7 cm~ 3 that we 
infer for BN, collisional broadening of recombination 
lines can be significant. Collisional broadening is most 
important for high-?? transitions because the diameter 
of a hydrogen atom is proportional to r? 2 - an atom 
in the n=53 level is 0.3 /im across! The collision- 
broad ened linewidths are given by the approximate for- 
mula (jBrocklehurst fc Seatonll 1972ft 

A?; c /2 = 4.7 e/v (??/100) 4 - 4 (10 4 /T e ) al r? e . (3) 

For n e = 5 x 10 7 cm" 3 and T e = 8000 K, Av c ~ 
200 km s" 1 for H53a, 30 km s" 1 for H41a, and 3 km s _1 
for H30a. These widths must be convolved with the ex- 
pected ~ 20 km s _1 wide thermal profiles, plus any ad- 
ditional Doppler broadening from turbulent motions. 

Since the observed width of the H53a line is much less 
than 200 km s _1 , most of the H53o; emission that is de- 
tected must originate more than ~ 12 AU from the star, 
in the lower density halo where n e < 10 7 cm -3 . The line 
is slightly blueshifted relative to the other transitions, as 
expected if this gas is expanding. The H30a line is little 
affected by pressure broadening or by continuum opacity, 
so its central velocity, 23.2 km s _1 , is the best indicator 
of the Vlsr of BN. We attribute the anomalously narrow 
width observed for the H41/42a transition to poor signal 
to noise; more sensitive spectra of all these recombina- 
tion lines would be useful in further constraining models 
of the electron density in BN. 

We do not detect the 231.995 GHz helium recombi- 
nation line toward BN. The absence of this line is con- 
sistent with the identification of BN as a B-star. In the 
30 km s _1 resolution BN spectrum in Figure[IJ the 2a up- 
per limit on the intensity of the Hc30a line is ~ 12 mjy, 
while the intensity of the H30a line is 270 mjy; thus 
r(Hc30a)/r(H30a) < 0.045. If helium were ionized over 
the full volume of the HII region, the opacity ratio would 
be comparable to the He/H abundance ratio, ~ 0.1. 
Thus we can say that the He + zone, if it exists, fills less 
than half the volume of the HII region. This places an up- 
per limit of roughly 0.05 on the ratio of heli um- ionizing to 
hydro gen- ionizing photons from the star (jMezger et al.1 
119741 Fig. 2), consistent with the ratio ex pected for B- 
stars , but not for 09 or more massive stars (IMezger et al.1 
IWl^Nacca et al.|[l996h . 

4.2. Srcl 

According to the data in Figure [31 the flux density 
of Srcl scales approximately as v 2 from 43 to 229 GHz, 
suggesting that over this frequency range the source is 
an optically thick black body with constant angular size. 
Even if the absolute calibrations arc wrong, the absence 
of a detectable H30a line rules out the possibility that the 
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Table 3 

Recombination Lino Parameters for BN 



line 


frcq 
(GHz) 


o a 
'-'peak 

(mJy) 


Vlsr 
(km s — 1 ) 


FWHM 
(km s — 1 ) 


H53a b 
H42o 
H41<* 
H30a 


42.952 
85.688 
92.034 
231.901 


10.4 ±1.1 
47±5 C 
291 ± 10 


20.1 ± 2.1 

24.2 ± l.l c 
23.2 ±0.5 


39.0 ±4.9 
20.0 ± 2. 7 C 
31.8 ± 1.3 



a Peak flux density of the line above the continuum level. 
b IRodriguez et al.ll2009l . 

c fit to average of the H41a and H42o line profiles. 



emission is optically thin p + /c~ bremsstrahlung, unless 
the H30a line is unexpectedly broad (> 300 km s" 1 ). For 
a linewidth of 30 km s~ 4 , comparable to the velocity span 
of the Srcl SiO masers, the line to continuum opacity 
ratio tl/tq ~ 3 from equations (1) and (2), and the line 
should easily have been detected if tc < 2. 

Our data cannot easily be reconciled with the 348 GHz 
flux density of 320 ± 48 mJy measured with the SMA 
(|Beuther et al.ll2~004l ). which falls below plausible extrap- 
olations from the lower frequency points. Such a low 
348 GHz flux density would imply that the free- free emis- 
sion from Srcl is somewhat optically thin at 229 GHz, in 
which case we should have detected the H30a line. Pos- 
sibly the 348 GHz Srcl measurement was corrupted by 
negative sidelobes from the bright Hot Core clumps 1-2" 
east of Srcl, or possibly Srcl is time-variable (cf. Sec- 
tion • The 690 GHz flux density me asured with the 
SMA, 6.7 ± 3.2 Jy (jBeuther et al.l 120061 ). lies above the 
v 2 extrapolation from lower frequencies, but this is easily 
explained if there is thermal emission from dust within 
the synthesized beam. 

4.2.1. An optically thick hypercompact HII region? 

If Srcl is a conventional hypercompact HII region, 
it must have an extraordinarily high emission measure, 
> 3 x 10 11 pc cm -6 , in order to remain optically thick 
at 229 GHz. We are able reproduce the observed contin- 
uum flux densities with the simple spherically symmetric 
model described in section 4.1.1 if we assume that the 
ionized region has radius r c = 7.5 AU, electron tem- 
perature T e = 8000 K, and uniform electron density 
n e = 1.5 x 10 8 cm -3 . The resulting continuum spec- 
trum is shown by the solid curve in Figure |3j the dashed 
curve shows the recombination line intensities predicted 
by the model. The HII region is constrained to have 
a sharp outer edge - any halo of lower density plasma 
would lead to excess emission at cm wavelengths. Col- 
lisional broadening of the H30a line is predicted to be 
of order 10 km s _1 , not enough to smear out the line. 
The excitation parameter U ~ 10 pc cm~ 2 computed 
from the radius and electron density corresponds to an 
ionizing flux of 4 x 10 46 photons s _1 , consistent with a 
B0 to Bl main sequence star wit h luminosity ~ 10 4 Lq 
(jPanagial 11975 ISmith "eTaLl 120021 ) . Since the continuum 
emission is optically thick, this should be interpreted as 
a lower limit to the ionizing flux. 

The source diameter of 15 AU used in the model is 
dictated by the assumed electron temperature of 8000 K, 
the measured flux densities, and the distance to Orion. A 



diameter of 15 AU is, however, inconsistent with the ob- 
served size of Srcl, which is roughly 0.2", or 80 AU, along 
its major axis. This suggests that the emission region is 
clumpy or filamentary so that it does not fill the synthe- 
sized beam. But it is implausible that all these clumps or 
filaments are perfectly sharp-edged, such that every one 
of them is optically thick even at 229 GHz. Generally, 
one expects a distribution of clump s to produce emission 
with a shallower spectral slope ([Ignace fe Churchwelll 
[200l . 

One alternative is that the electron temperature is less 
than 8000 K. The highest resolution measurements of 
Srcl that are available, 43 GHz VLA observations with 
a 34 mas synthesized beam, imply tha t the source has a 
brightness temperature of only 1500 K dReid et al.ll2007D . 
This is much cooler than is possible in an HII region, 
however; hydrogen is almost entirely n eutral in dense 
gas at T< 4000 K (|Reid fc Mentedll997l ). 

4.2.2. H~ free-free emission from a massive disk? 

The absence of hydrogen recombination lines, 1500 K 
apparent brightness temperature, and spectral index of 
~2 all are difficult to explain if Srcl is a hypercom- 
pact HII region. The data more nat urally fit an alter - 
native hypothesis, first considered bv lReid et al.l (|2007| ). 
that the emission from Srcl originates in a T< 4500 K 
disk, similar to the radio photospheres of Mira vari- 
ables (|Reid fe Mentenlll997f ). The infrared spectrum of 
Srcl, seen via light reflected off the neighbori ng nebulos- 
ity, su ggests such an interpretation as well - iTesti et ail 
(|2010T l find that the spectrum is characteristic of a low- 
gravity photosphere with effective temperature 3500- 
4500 K. The close association of Srcl with SiO masers, 
also characteristic of Mira variables, further strengthens 
the case for this h ypothesis. 

As discussed bv lReid fc Mentenl ()1997f ). free- free emis- 
sion at temperatures of about 1500 K in such a photo- 
sphere arises via the H~ opacity, the scattering of elec- 
trons by neutral hydrogen atoms or molecules. The elec- 
trons are produced by collisional ionization of Na, K, 
and other metals. The H~ opacity is roughly 10 3 times 
smaller than the p + /e~ opacity, so high densities (> 10 11 
cm~ 3 ) typically are required to obtain optically thick 
e mission. 

iReid et al.1 (|2007f ) were able to fit the 43 GHz VLA ob- 
servations of Srcl with models of H~ emission from an 
80 AU diameter, 3 M© disk. However, they found that 
such a large disk cannot be maintained at T > 1500 K 
if it is heated solely by a central star - a stellar lumi- 
nosity > 10 5 L/0, greater than the luminosity of the 
entire Orion-KL region, would be required. They pro- 
posed that accre t ion p rocesses might heat the disk lo- 
cally. ITesti et aTl (|2010D also find that the infrared spec- 
trum of Srcl can be explained by emission from a disk 
around a ~ 10 M© protostar that is accreting at a 
few x 10 -3 Mq yr _1 . However, their models show that 
the disk temperature falls below 1500 K for r> 10 AU, so 
cannot easily explain the H~ free-free emission at larger 
radii. 

In summary, the H~ hypothesis explains the observa- 
tional results in a natural way. The principal barrier 
to its acceptance is uncertainty over the mechanism by 
which such a large disk can be maintained at sufficiently 
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Figure 5. SrcI/BN flux density ratio vs. year, at 43 GHz (filled 
squares) and at 86 GHz (circles), derived from the data in Table[2] 
This ratio should be immune to errors in the absolute flux scale, 
since the Srcl and BN arc just 10" apart on the sky and are ob- 
served simultaneously. These data provide possible evidence that 
Srcl has brightened relative to BN since 1995. 

high temperature. Whether the disk could be heated by 
an anomalously high accretion rate, a merging binary, 
magnetic effects, or some other mechanism remains a key 
question. 

4.2.3. Is Srcl brightening? 

If accretion is responsible for much of SrcFs luminosity, 
then one might expect to see c hanges in its flux den sity on 
timescales of years. Although |(Toddi et all (pOlTf ) found 
no evidence of flux variations at 43 GHz over 4 epochs 
from 2000 to 2009, we are struck by the discrepancy in 
the 86 GHz fluxes mea sured at BIMA in 1995 (34±5 mJy, 
iPlambeck etTaTlll995l) and at CARM A in 2009 (50 ± 5 
mJv. lFriedel k Widicus Weavenl2^TlT ). 

Generally, uncertainties in absolute flux densities are 
dominated by systematic effects - e.g., variations in tele- 
scope gains as a function of elevation, or uncertainties 
in the brightness temperature models of planets used as 
primary flux standards. The ratio of the flux densities of 
Srcl and BN should be immune to such effects because 
the two sources are just 10" apart and can be observed 
simultaneously. In Figure [5] we plot this ratio as a func- 
tion of time for the 43 and 86 GHz data in Table [2j Note 
that the error bars in this figure were calculated from 
the thermal noise in the original maps, not from the ab- 
solute flux uncertainties listed in Table [2j From these 7 
independent datasets it appears that the ratio increased 
from about 0.4 in 1995 to 0.5 in 2009. 

This result should be interpreted cautiously. Al- 
though the flux density ratio is unaffected by the ab- 
solute calibration scale, it is sensitive to the sampling 
of visibility data in the (u, v) plane. It is also unclear 
whether Srcl or BN, or both, are varying; periodic vari- 
ations in the infra r ed lum inosity of BN were reported by 
iHillenbrand et aD (|2001h . Future careful monitoring of 
the flux densities of both Srcl and BN with the VLA 
and ALMA will be required to find conclusive evidence 
for source brightness variations. 

4.2.4. Vibrationally excited H% O 

The bright emission line in the spectrum of Srcl in 
Figure 1 is the H 2 v 2 =l, 5(5,0)-6(4,3) transition of H 2 
at 232.687 GHz, 3450 K above the ground state. This 



line also was recognized by iHirota et al.l (|2012[ ) in the 
Orion Band 6 ALMA Science Verification data. With 
0.6 km s _1 velocity resolution, the ALMA data show 
that the H 2 line is double-peaked, like the 22 GHz H 2 
masers or the 43 and 86 GHz SiO masers toward Srcl. 

In the 1.5" resolution ALMA maps the H 2 line 
is blended with a HCOOCH 3 transition at 232.684 
GHz, whereas in the 0.25" resolution CARMA map the 
HCOOCH3 emission is resolved out. The H 2 emission 
is unresolved by the CARMA beam, so it is clear that the 
line originates from the compact SiO maser zone close to 
Srcl, rather than the extende d 2" x 0.5" strip wh ere 22 
GHz H 2 masers are found (|Gaume et al.lll998D . The 
SiO v=2 masers also ori ginate in energy lev els 3500 K 
above the ground state; iGoddi et aTl ( 20091) find that 
these masers can be excited in gas with kinetic temper- 
ature > 2000 K. 

The same transition of vibrationally excited H 2 was 
det ected toward the evolved stars VY CMa and W Hya 
bv iMenten k Melnickl (|1989fl . Like Srcl, both of these 
stars have SiO and 22 GHz H 2 masers. 

5. SUMMARY 

CARMA 229 GHz maps of Orion-KL with up to 0.14" 
angular resolution clearly distinguish the ionized en- 
velopes of Srcl and BN from the surrounding molecular 
cloud. The principal results arc as follows: 

1 . The integrated flux densities of Srcl and BN at 229 
GHz are 310 and 240 mJy respectively, with a prob- 
able uncertainty of ±15%. The BN flux density we 
measure with 0.14" resoluti on is a factor of two 

reater than that derived bv iGalvan-Madrid et al.l 
20121 ) from ALMA Science Verification data with a 
1.3" x 0.6" synthesized beam. We suspect that the 
flux density of BN is depressed in the ALMA map 
because of negative sidelobes from poorly sampled 
extended emission. These large scale structures are 
filtered out more effectively in the higher resolution 
CARMA data. 

2. Srcl's flux density is proportional to v 2 from 43 to 
229 GHz - it appears to be an optically thick black 
body. By comparison, BN's flux density scales as 
v 1 ' 3 up to about 100 GHz, then flattens; its free-free 
continuum is becoming optically thin at 229 GHz. 

3. Srcl is elongated at PA 140 degrees in the CARMA 
images. Its deconvolved size is similar at 229 and 
43 GHz - further evidence that it is optically thick. 

4. The H30a recombination line is not detected to- 
ward Srcl, though it is bright toward BN. The ab- 
sence of this recombination line is a third piece of 
evidence that Srcl is optically thick, or that it is 
not an HII region at all. 

5. The v 2 =l, 5(5,0)-6(4,3) transition of H 2 at 
232.687 GHz, 3450 K above the ground state, is 
prominent toward Srcl. The line emission origi- 
nates close to the central star, not in the larger 
zone where 22 GHz H 2 masers are found. 

6. While our observational results do not rule out the 
possibility that Srcl is an exceptionally dense hy- 
percompact HII region, they are more easily ex- 
plained by models of free-free emission via the 
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H opacity in a T< 4500 K disk, as modeled by 
IReid etaLl (12001 . 

7. A model of a spherically symmetric hypercompact 
HII region with a diameter of about 15 AU and 
electron density n e ~5x 10 7 cm -3 matches both 
the continuum and the recombination line intensi- 
ties toward BN. 

8. The excitation parameter U ~ 5 pc enr 2 derived 
for BN implies that its Lyman continuum flux is 
5 x 10 45 photons s _1 , consistent with that expected 
from a zero age main sequence B star with a total 
luminosity of (0.5 - 1) x 10 4 L Q . 

9. The H53a recombination line in BN is predicted 
to have a collisionally broadened lincwidth of 
200 km s -1 . Presumably the H 53a emission ob- 
served bv lRodriguez et all (|2009l) . with a linewidth 
of about 40 km s , originates primarily from lower 
density ionized gas farther from the central star. 

10. The V LS r of BN is 23.2 ± 0.5 km s" 1 , based on 
the central velocity of the 231.9 GHz H30a line, 
which is little affected by pressure broadening or 
continuum opacity. 

11. Published flux densities at both 43 and 86 GHz pro- 
vide tentative evidence that Srcl brightened with 
respect to BN between 1994 and 2009. 

Future ALMA observations of Orion-KL at submil- 
limeter wavelengths will provide a definitive test of the 
emission mechanism in Srcl. If Srcl is an HII region, it 
must become optically thin at sufficiently high frequen- 
cies. If the source were still optically thick at 650 GHz, 
for example, an O-star with a Lyman continuum flux 
> 10 49 photons s _1 would be required to maintain its 
ionization, but the luminosity of such a star, > 10 5 Lq 
(|Panagiall973[ ). would exceed the measured luminosity of 
the entire Orion-KL nebula. Then, at frequencies where 
the continuum is optically thin, hydrogen recombination 
lines must be detectable if hydrogen is ionized. If these 
recombination lines arc detected, the free-free hypothesis 
is confirmed. If not, or if the continuum flux density con- 
tinues to increase following a v 2 law, the emission must 
be attributable to the H~ mechanism. 

Submillimctcr observations of BN, as well as more sen- 
sitive observations of its cm and mm wavelength recom- 
bination lines, will be useful in constraining models of 
the hypercompact HII region around this star, and in 
tracking small changes in the SrcI/BN flux density ra- 
tio that could be an indication of time variable accretion 
onto these objects. 
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tional Science Foundation. Ongoing CARMA develop- 
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ence Foundation under a cooperative agreement, and by 
the CARMA partner universities. 
Facilities: CARMA. 

REFERENCES 

Allen, D. A., & Burton, M. G. 1993, Nature, 363, 54 

Bally, J., Cunningham, N. J., Moeckel, N., Burton, M. G., Smith, 

N., Frank, A., & Nordlund, A. 2011, ApJ, 727, 113 
Beuther, H., et al. 2004, ApJ, 616, L31 
— . 2006, ApJ, 636, 323 

Brockletmrst, M., & Seaton. M. J. 1972. MNRAS, 157, 179 
Chandler, C. J., & Wood, D. O. S. 1997, MNRAS, 287, 445 
De Buizer, J. M., Morris, M. R., Becklin, E. E., Zinnecker, H., 

Herter, T. L., Adams, J. D., Shuping, R. Y., & Vacca, W. D. 

2012, ApJ, 749, L23 
Doi, T., O'Dell, C. R., & Hartigan, P. 2002, AJ, 124, 445 
Eisner, J. A., Plambcck, R. L., Carpenter, J. M., Corder, S. A., 

Qi, C, & Wilner, D. 2008, ApJ, 683, 304 
Felli, M., Churchwell, E., Wilson, T. L., & Taylor, G. B. 1993a, 

A&AS, 98, 137 

Felli, M., Taylor, G. B., Catarzi, M., Churchwell, E., & Kurtz, S. 

1993b, A&AS, 101, 127 
Forbrich, J., Menten, K. M., & Reid, M. J. 2008, A&A, 477, 267 
Friedel, D. N., & Widicus Weaver, S. L. 2011, ApJ, 742, 64 
Galvan-Madrid, R., Goddi, C, & Rodriguez, L. F. 2012, A&A, 

547, L3 

Garay, G., Moran, J. M., & Reid, M. J. 1987, ApJ, 314, 535 
Gaume, R. A., Wilson, T. L., Vrba, F. J., Johnston, K. J., & 

Schmid-Burgk, J. 1998, ApJ, 493, 940 
Goddi, C, Greenhill, L. J., Chandler, C. J., Humphreys, 

E. M. L., Matthews, L. D., & Gray, M. D. 2009, ApJ, 698, 1165 
Goddi, C, Humphreys, E. M. L., Greenhill, L. J., Chandler, 

C. J., & Matthews, L. D. 2011, ApJ, 728, 15 
Gomez, L., Rodriguez, L. F., Loinard, L., Lizano, S., Allen, C, 

Poveda, A., & Menten, K. M. 2008, ApJ, 685, 333 
Hillenbrand, L. A., Carpenter, J. M., & Skrutskie, M. F. 2001, 

ApJ, 547, L53 

Hirota, T., Kim, M. K., & Honma, M. 2012, ApJ, 757, LI 
Ignace, R., & Churchwell, E. 2004, ApJ, 610, 351 
Kim, M. K., et al. 2008, PASJ, 60, 991 

Matthews, L. D., Greenhill, L. J., Goddi, C, Chandler, C. J., 

Humphreys, E. M. L., & Kunz, M. W. 2010, ApJ, 708, 80 
Menten, K. M., & Melnick, G. J. 1989, ApJ, 341, L91 
Menten, K. M., & Reid, M. J. 1995, ApJ, 445, L157 
Menten, K. M., Reid, M. J., Forbrich, J., & Brunthaler, A. 2007, 
A&A, 474, 515 

Mezger, P. G., Smith, L. F., & Churchwell, E. 1974, A&A, 32, 269 
Morino, J. -I., Yamashita, T., Hasegawa, T., & Nakano, T. 1998, 

Nature, 393, 340 
Panagia, N. 1973, AJ, 78, 929 
Perez, L. M., et al. 2010, ApJ, 724, 493 
Plambeck, R. L., et al. 2009, ApJ, 704, L25 
Plambeck, R. L., Wright, M. C. H., Mundy, L. G., & Looney, 

L. W. 1995, ApJ, 455, L189+ 
Reid, M. J., & Menten, K. M. 1997, ApJ, 476, 327 
Reid, M. J., Menten, K. M., Greenhill, L. J., & Chandler, C. J. 

2007, ApJ, 664, 950 
Rodriguez, L. F., Poveda, A., Lizano, S., & Allen, C. 2005, ApJ, 

627, L65 

Rodriguez, L. F., Zapata, L. A., & Ho, P. T. P. 2009, ApJ, 692, 
162 

Rohlfs, K., & Wilson, T. L. 2000, Tools of Radio Astronomy (3rd 

edition; New York: Springer) 
Schaller, C, Schaerer, D., Mcynet, G, & Maeder, A. 1992, 

A&AS, 96, 269 

Scoville, N., Kleinmann, S. G., Hall, D. N. B., & Ridgway, S. T. 

1983, ApJ, 275, 201 
Smith, L. J., Norris, R. P. F., & Crowther, P. A. 2002, MNRAS, 

337, 1309 

Testi, L., Tan, J. C, & Palla, F. 2010, A&A, 522, A44 
Vacca, W. D., Garmany, C. D., & Shull, J. M. 1996, ApJ, 460, 914 
Wright, A. E., & Barlow, M. J. 1975, MNRAS, 170, 41 
Zapata. L. A., Loinard, L., Schmid-Burgk, J., Rodriguez, L. F., 
Ho, P. T. P., & Patel, N. A. 2011, ApJ, 726, L12 



